## Not run:
# pompExample(ou2)
#
# guess1 <- guess2 <- coef(ou2)
# guess1[c('x1.0','x2.0','alpha.2','alpha.3')] <- 0.5*guess1[c('x1.0','x2.0','alpha.2','alpha.3')]
# guess2[c('x1.0','x2.0','alpha.2','alpha.3')] <- 1.5*guess1[c('x1.0','x2.0','alpha.2','alpha.3')]
#
# m1 <- mif2(ou2,Nmif=100,start=guess1,Np=1000,
# cooling.type="hyperbolic",cooling.fraction.50=0.05,
# rw.sd=rw.sd(x1.0=ivp(0.5),x2.0=ivp(0.5),
# alpha.2=0.1,alpha.3=0.1))
#
# m2 <- mif2(ou2,Nmif=100,start=guess2,Np=1000,
# cooling.type="hyperbolic",cooling.fraction.50=0.05,
# rw.sd=rw.sd(x1.0=ivp(0.5),x2.0=ivp(0.5),
# alpha.2=0.1,alpha.3=0.1))
#
# plot(c(m1,m2))
#
# rbind(mle1=c(coef(m1),loglik=logLik(pfilter(m1,Np=1000))),
# mle2=c(coef(m2),loglik=logLik(pfilter(m1,Np=1000))),
# truth=c(coef(ou2),loglik=logLik(pfilter(m1,Np=1000))))
# ## End(Not run)
Run the code above in your browser using DataLab